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ABSTRACT 


In  the  Low  Frequency  (LF,  30-300  kHz)  band  of  the 
electromagnetic  spectrum,  one  finds  many  digital  commu- 
nication systems.   At  these  frequencies,  atmospheric 
noise,  predominantly  caused  by  radiation  from  lightning 
discharges,  is  the  dominant  noise  which  corrupts  the  sig- 
nal.  When  viewed  at  the  output  of  the  receiver's  band- 
limiting  filter,  the  noise  waveform  is  very  definitely 
non-gaussian.   Since  the  noise  is  non-gaussian,  it  would 
be  difficult  to  design  an  optimum  communication  receiver 
without  knowledge  of  specific  properties  of  the  noise. 

A  model  for  bandlimited  low  frequency  atmospheric 
noise,  which  accurately  fits  both  the  first  order 
probability  density  and  the  time  structure  of  observed 
noise  waveforms,  and  which  is  well  suited  for  digital  ^ 
computer  simulations,  has  been  developed  by  Feldman. 
In  the  development  of  this  model,  a  statistical  dependence 
between  the  noise  waveform  on  one  bandlimited  frequency 
channel  and  the  envelope  of  the  noise  on  a  disjoint,  but 
nearby,  frequency  channel  was  identified  and  incorporated 
into  the  model.   On  the  basis  of  this  property,  a  design 
for  a  two  channel  low  frequency  communication  receiver, 
which  makes  nearly  optimum  use  of  the  joint  channel 
information,  was  proposed. 
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This  work  has  addressed  itself  to  the  task  of 
simulating  the  joint  channel  receiver  and  testing  its 
performance  relative  to  a  single  channel  linear  corre- 
lation receiver.  Both  receivers  operate  in  the  presence 
of  low  frequency  atmospheric  noise,  as  simulated  by  the 
above  mentioned  noise  model. 

The  modem  was  assumed  to  be  coherent  phase  shift 
keying,  with  the  signal  set  consisting  of  ±   A  cos  2nf  t. 
The  criterion  for  performance  comparison  was  an  estimate 
of  probability  of  error  per  binary  bit  of  information 
received.   Our  results  show  that,  with  the  dual  channel 
receiver,  performance  improvements  of  more  than  an  order 
of  magnitude  can  be  expected  during  periods  of  severe 
non-gaussian  noise  activity.   Suggestions  for  further 
research  with  this  receiver  design  are  presented  in  the 
last  chapter  of  this  work. 
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Chapter  1 

INTRODUCTION 

In  the  Low  Frequency  (LP,  30-300  kHz)  band  of  the 
electromagnetic  spectrum,  one  finds  many  digital  commu- 
nication systems.   At  these  frequencies,  the  dominant 
noise  encountered  is  atmospheric  noise,  which  is 
predominantly  caused  by  radiation  from  lightning 
discharges.   Since  individual  discharges  tend  to 
dominate  the  noise  waveform  observed  at  the  output 
of  a  bandlimiting  filter,  we  must  assume  that  the  noise 
cannot  be  statistically  described  as  a  gaussian  noise 
process.   Clearly  then,  one  cannot  hope  to  build  an 
optimum  communication  receiver  without  some  knowledge 
of  the  specific  non-gaussian  characteristics  of  the  low 
frequency  atmospheric  noise. 

The  research  reported  in  this  work  will  demonstrate 
one  particular  realization  of  a  near  optimum  binary 
digital  communication  receiver  which  is  based  upon  the 
observed  characteristics  of  bandlimited  atmospheric 
noise.   In  this  chapter,  we  will  discuss  the  actual 
atmospheric  noise  and  a  model  for  that  noise  which 
can  be  used  in  digital  computer  simulations. 

In  Chapter  2,  we  present  a  receiver  design  which 
uses  a  second  disjoint  frequency  channel  to  provide  more 
information  about  the  noise  which  corrupts  the  signal. 
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We  will  show  that  this  receiver  makes  nearly  optimum  use 

of  the  joint  channel  information.   We  will  also  compare 
this  receiver  to  a  single  channel  linear  correlation 
receiver. 

Chapter  3  is  devoted  to  an  explanation  of  the 
calculations  necessary  to  relate  the  analog  system  to 
the  digital  simulations  of  the  proposed  receiver.   Our 
criterion  for  evaluation  of  receiver  performance  will 
simply  be  an  estimate  of  the  probability  of  error  per 
bit  of  binary  information  received. 

In  Chapter  4,  we  will  present  the  results  of  our 
simulations  of  both  the  single  channel  and  joint  channel 
receiver  models,  in  the  presence  of  the  simulated  low 
frequency  atmospheric  noise.   We  will  give  bounds  on 
the  accuracy  of  our  results  and  explain  some  of  the 
problems  encountered  in  carrying  out  the  simulations. 

Finally,  Chapter  5  will  give  a  summary  of  this 
research,  the  conclusions  reached,  and  some  suggestions 
for  further  research, 

1.1  Low  Frequency  Atmospheric  Noise 

There  have  been  many  observations  of  bandlimited 
low  frequency  atmospheric  noise,  which  have  led  to 
descriptions  of  the  noise  characteristics  (for  example, 
see  1,2),  but  all  note  that  the  noise  waveform  is  an 
extremely  non-gaussian  process.   At  individual  times, 
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the  noise  waveform  Is  dominated  by  peaks  of  high 
instantaneous  energy,  with  several  peaks  occasionally 
occurring  rapidly  enough  to  yield  a  seemingly  continuous 
burst  of  high  energy. 

Each  high  energy  peak  can  be  related  to  the 
radiation  released  from  a  lightning  discharge.   There 
are  two  basic  strokes  in  each  discharge:   the  leader 
stroke  with  an  average  length  of  about  1  msec  and 
centered  in  frequency  about  30  kHz,  and  the  return 
stroke  with  a  length  of  about  100-200  usee  and  centered 
about  3-8  kHz.   In  many  instances,  there  are  complicated 
repititions  of  these  two  basic  strokes  which  yield  the 
continuous  noise  burst  behavior  noted  above.   Between 
high  energy  bursts,  the  noise  waveform  can  be  accurately 
characterized  as  a  low  level  gaussian  waveform. 

Feldman  reported  that  the  following  observed 
properties  of  the  bandlimited  noise  waveform  were  useful 
in  developing  a  model  for  this  noise:   First,  at  low 
frequencies,  the  noise  is  highly  dependent  upon  local 
weather  conditions.   Second,  above  a  certain  threshold, 
the  first  order  joint  statistics  of  the  noise  waveform 
in  one  frequency  channel  and  the  noise  envelope  in  a 
second  disjoint,  but  nearby,  frequency  channel  are 
statistically  dependent.   Third,  the  bandlimited  waveform 
has  an  average  power  spectral  density  which  is  proportional 
to  the  square  of  the  narrow-band  filter's  frequency 
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response.   And  fourth,  for  bandwidths  of  1  kHz  or  larger, 
the  process  can  be  considered  uncorrelated  beyond  1  msec. 
In  the  next  section,  we  will  present  the  noise  model  which 
Feldman  subsequently  developed. 

In  his  analysis  of  atmospheric  noise,  Feldman  noted 
that  one  of  four  basic  weather/noise  conditions  is 
sufficient  to  describe  the  noise  at  any  particular  time. 
These  noise  conditions  are  quiet,  tropical,  frontal,  and 
quiet  night.   Quiet  corresponds  to  a  stable  air  mass, 
with  little  thunderstorm  activity,  in  the  vicinity  of  the 
receiver.   Tropical  conditions  exist  when  warm,  moist, 
unstable  air,  with  much  thunderstorm  activity,  is  present 
near  the  receiver.   Frontal  corresponds  to  the  passage 
of  a  front  in  the  vicinity  of  the  receiver.   Quiet 
night  represents  a  transitional  period  between  quiet  and 
tropical  conditions.   The  order  of  severity  varies  from 
nearly  gaussian  to  nearly  continuous  high  energy  burst 
behavior  for  quiet,  quiet  night,  tropical,  and  frontal 
conditions,  respectively. 

1.2  Noise  Model 

The  model  for  bandlimited  low  frequency  atmospheric 
noise,  y(t),  which  we  will  consider  in  this  work  is  the 
following: 

y(t)  =  ni(t)  +  a(x(t))n2(t),        (1.1) 
where  n^(t)  and  ^(t)  are  statistically  independent 
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gaussian  processes  and  n-j_(t)  represents  the  low  level 
gausslan  background  process.   The  quantity  a(x(t))n2(t) 
represents  the  high  energy  bursts  of  noise  due  to 
lightning  discharges,  and  x(t)  is  a  two-state  Markov 
process  with  a  time-varying  rate  parameter.   The  two 
states  of  x(t)  turn  on  and  off  the  high  energy  noise 
component.   Whenever  x(t)  is  in  the  "off"  state,  a(.) 
is  zero,  and  whenever  x(t)  is  in  the  "on"  state,  a(.) 
becomes  a  random  number,  which  is  constant  over  the 
time  interval  that  x(t)  remains  in  the  "on"  state. 
Each  time  x(t)  changes  from  "off"  to  "on,"  a  new  random 
value  for  a(.)  is  generated,  and  a(.)  remains  constant 
until  x(t)  turns  it  off. 

This  model  incorporates  the  observed  dependence 
between  the  noise  waveform  on  one  channel  and  the 
envelope  of  the  noise  on  a  second  disjoint,  but  nearby, 
frequency  channel  by  assuming  that  a(x(t))  remains 
independent  of  frequency.   Thus,  it  follows  that  the 
model  for  the  envelope  of  the  bandlimited  noise,  v(t), 
on  a  second  frequency  channel  is 

v(t)  =  Vl(t)  +  a(x(t))v2(t),        (1.2) 
where  v-^(t)  and  v2(t)  are  statistically  independent 
Rayleigh  processes. 

It  must  be  noted  that  this  noise  model  cannot 
exactly  specify  the  short-time  structure  of  the  actual 
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atmospherlc  noise*  However,  since  LF  receivers  use  sharp 

band limiting  filters,  with  long  observation  times,  this 

noise  model  is  sufficient  because  it  accurately  simulates 

the  average  behavior  of  the  noise  over  these  relatively 

long  tine  intervals.   The  model  is  based  on  the  first 

order  probability  density  of  the  bandlimited  noise.   In 

addition,  properties  of  a(x(t))  are  such  that  the  model 

accurately  simulates  the  actual  long-time  time  structure 

of  the  noise. 

Each  of  the  four  noise  conditions  (quiet,  quiet  night, 

tropical,  and  frontal)  can  be  simulated  with  the  noise 

model  simply  by  varying  the  parameters  used  in  the 

generating  mechanism.   In  this  research,  all  four  noise 

conditions  were  simulated  extensively  using  a  canonic 

noise  generator  and  table  of  associated  parameters 

given  in  Figure  3-H  and  Table  3-1 >  respectively,  of 

reference  (1).   Since  these  two  items  are  of  vital 

importance  to  this  research,  they  have  been  included 

as  Figure  1.1  and  Table  1.1  of  this  report.   The  results 

of  the  simulations,  conducted  in  this  research  project, 

are  given  in  Chapter  4  of  this  work. 
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Chapter  2 
RECEIVER 

In  Chapter  1,  we  presented  some  properties  of  the 
bandlimited  low  frequency  atmospheric  noise,  the  noise 
which  corrupts  signals  in  low  frequency  radio  communication 
systems.   We  also  presented  a  model,  suitable  for  simulation 
on  a  digital  computer,  which  accurately  represents  that 
noise.   The  low  frequency  atmospheric  noise,  when  viewed 
at  the  output  of  the  receiver's  bandlimiting  filter,  has 
been  observed  to  be  extremely  non-gaussian.  However,  we 
presented  several  properties  of  the  noise  which  can  be 
exploited  in  the  design  of  a  nearly  optimum  receiver. 

In  this  chapter,  we  will  present  a  model  for  a 
phase  coherent  linear  correlation  receiver.  We  will 
show  that  this  is  the  optimum  receiver  for  signals 
corrupted  by  white  gaussian  noise,  but  it  is  not  neces- 
sarily optimum  when  the  noise  is  non-gaussian.   We  will 
also  present  a  receiver  design,  based  on  the  specific 
observed  properties  of  bandlimited  low  frequency 
atmospheric  noise,  which  should  yield  nearly  optimum 
performance.  However,  implementation  of  the  new  design 
is  approximately  twice  as  complicated  as  implementation 
of  the  linear  correlation  receiver. 
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2.1  Modulation-Demodulation  Scheme 

Since  the  basic  purpose  of  this  research  is  to 
compare  the  performance  of  a  new  receiver  design  with 
that  of  a  standard  receiver  design,  we  have  chosen  a 
computationally  simple  modulation-demodulation  scheme. 
We  postulate,  since  we  are  comparing  both  receivers 
under  the  same  set  of  conditions,  that  our  results  in 
this  somewhat  idealistic  case  will  still  be  representative 
of  the  results  we  would  expect  in  a  more  realistic  case. 

In  the  simplest  case,  the  receiver  must  decide 
whether  or  not  a  signal  is  present.   In  a  slightly  more 
complicated,  and  more  realistic,  case,  the  receiver  must 
decide  between  two  signals.   We  have  thus  chosen  our  sig- 
nal set  to  contain  two  signals,  each  containing  the  same 
amount  of  energy.   This  signalling  scheme,  in  which  one 

signal  is  the  negative  of  the  other,  is  known  as  anti- 

3 
podal  signalling. 

Since  only  a  limited  amount  of  data  was  anallzed  in 

the  development  of  the  noise  model,  we  have  parameters 

available  to  model  the  bandlimited  noise  waveform  only 

under  certain  frequency  and  bandwidth  constraints. 

Keeping  these  constraints  in  mind,  we  must  satisfy  the 

not  unrealistic  condition  that  all  radiated  signal 

energy  be  confined  to  the  frequency  band  centered  at 

65  kHz,  with  a  bandwidth  of  1  kHz.1  (see  Table  1.1) 
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This  means  that  in  the  baseband,  the  lowpass  region  of 
the  frequency  spectrum,  our  signals  must  contain  no 
energy  in  any  frequency  component  greater  than  1  kHz, 
In  order  to  satisfy  this  constraint  in  a  simple  manner, 
we  have  chosen  our  baseband  signal  set  to  consist  of 

sb(t)  =  +  A,  (2.1) 

and  thus,  since  A  is  a  constant,  the  frequency  spectrum 
of  each  of  our  baseband  signals  consists  only  of  an 
impulse  at  the  frequency  origin. 

However,  we  are  restricted  to  a  frequency  band 
centered  not  at  the  origin  but  at  65  kHz  for  transmission. 
Therefore,  we  must  modulate  our  signals  up  to  this  pass- 
band.   To  do  this,  we  simply  allow  our  signals  to  ride 
on  a  cosine  carrier  of  frequency  65  kHz;  and  thus  our 
passband  signal  set  becomes 

s(t)   =  +  A  cos  2Ttfrtt,  (2.2) 

where  f     =  65  kHz. 

At  this  point,  we  will  separate  the  members  of  our 

transmitted  signal  set,  thus: 

s0(t)  =  A  cos  2rrf  t,  (2.3a) 

c 

which  represents  a  binary  zero  being  sent,  and 

s1(t)  =  -A  cos  277fct,  (2.3b) 

which  represents  a  binary  one  being  sent.   From  the  above, 
it  is  clear  that 

sj(t)  =  -sQ(t)  =  A  cos  (2iTfct  +  0),     (2.3c) 


-19- 

where  0   =  TT,   Thus,  we  see  that  our  signal  set  consists 

of  two  phase  shift  keyed  (PSK)  signals. 

In  general,  this  is  not  a  practical  way  of  signalling 
by  radio  communication  systems  due  to  the  high  voltage 
switching  transients  which  occur  when  reversing  phase. 
However,  for  signals  corrupted  by  gaussian  noise,  the 
optimum  PSK  demodulator  performs  3  db  better  than  the 
optimum  frequency  shift  keyed  (FSK)  demodulator,  where 
binary  FSK  signal  sets  consist  of 

sQ(t)   =  A  cos   2TTf0t  (2.4a) 

and 

s1(t)  =  A  cos  2fif1t,  f1   4   f0.       (2,4b) 
Furthermore,  low  frequency  signalling  schemes  which 
approach  PSK  in  performance  have  been  developed,  and  it 
is  felt  that  analysis  of  a  PSK  model  will  give  an 
Indication  of  performance  to  be  expected  in  these  new 
schemes. 

If  there  were  no  noise  present,  the  optimum  demodu- 
lator would  filter  out  all  signals  but  those  in  the 
passband,  monitor  the  received  signal,  compare  it  with 
a  reference  signal  (which  is  an  exact  replica  of  either 
sQ(t)  or  Sj(t)),  and  decide  which  signal  was  sent. 
However,  there  invariably  is  noise  present,  which  in  the 
simplest  form  is  additive  white  gaussian  noise.   Other 
forms  of  noise,  such  as  a  random  attenuation  of  the 
signal  or  a  random  phase  shift,  may  exist.   We  will, 
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however,  restrict  ourselves  to  the  additive  noise  case. 
The  reasons  for  this  will  be  given  in  the  next  section, 
where  we  explain  the  linear  correlation  receiver. 

2.2  Linear  Correlation  Receiver 

In  almost  any  text  on  communication  theory,  one 
can  find  a  description  of  a  linear  correlation  receiver 
(see,  for  example,  3»^*7).   In  Figure  2.1,  we  show  a 
linear  correlation  receiver  which  decides  which  of  two 
arbitrary  signals  has  been  received.   In  Figure  2.2,  we 
show  how  that  receiver  can  be  simplified  for  the  case  of 
antipodal  signals,  the  case  we  are  considering. 

Let  us  now  assume  that  the  received  signal,  r(t), 
consists  of  the  following  (see  Figure  2.2): 
r(t)  =  s(t)  +  n(t) 

=  +  A  cos  2-rrf^t  +  n(t),         (2.5) 
where  n(t)  is  the  additive  white  gaussian  noise.   We 
will  assume  that  n(t)  is  a  stationary,  zero  mean  pro- 
cess, with  autocorrelation  function, 

Rn(t1-t2)  =  ^0  %{t1't2),  (2.6) 

where  h ( . )  is  the  Dirac  delta  function.  From  this,  it 

is  clear  that 

rl(t)   =  ±  A2   cos2(2TTfct)   +  A  n(t)    cos   2lTfct 

=  +  A2    (1   +   cos  4T7f   t)   +  A  n(t)    cos   2irf„t,    (2.7) 
"  2~  c 

and 


-21- 


r(i) 


*& 


sAi) 


^ 


$- 


5.  ft) 


fc=r 


Figure  2,1  General  Linear  Correlation  Receiver 
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Figure  2.2  Linear  Correlation  Receiver- Antipodal  Signals 
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r?  =   1    f    f+  A2    (1  +  cos  ^-?rf„t)  +  A  n(t)  cos   2irf   t     dt 
TJOL"F                            c  CJ 

T 
=  +  A2  +   1    f      A  n(t)    cos   2TTf„t  dt,  (2.8) 

~  T      t-^o  c 

where  we  have  assumed  that  T,  the  bit  length  or  length 

of  time  in  which  the  receiver  must  decide  which  signal 

was  sent,  is  equal  to  an  integral  number  of  periods  of 

the  cosine  wave.   Actually,  this  is  not  a  restrictive 

assumption  because  the  bit  length  is  typically  much 

longer  than  the  length  of  one  period  of  the  cosine  wave. 

It  is  clear,  since  n(t)  has  been  defined  to  be  a 

gaussian  random  process,  and  since  r2  ^s  ^he  result  of 

a  linear  filtering  operation  on  n(t),  that  r2  is  a 

gaussian  random  variable.-'  Thus,  we  can  completely 

specify  the  probability  density  of  r2  if  we  can  find 

the  first  and  second  moments  of  r2. 

Taking  the  expected  value  of  equation  2.8,  we  get 

T 
E(r2)  =  +  a£  +   1.  E(f     A  n(t)    cos  2frfct  dt) 
"  2  T  0 

T 
■  ±  iL  +  i  f      A  E(n(t))    cos   2rTf   t  dt 
2  T^O  c 

-  +  •  (2.9) 


Taking  the  second  moment,  we  obtain 

T 
E(r?2)  =  (+  a£)2  +  a£  E(  C     2  A  n(t)  cos  2nf  t  dt) 
2      2T    0 
T  T 
+  L_  E(/  /  A2  n(t-)  n(t2)  cos  27tfct1 
T2    0   0 

cos  2rrfct2  dt.,  dt2) 


-  (±  a!)2 
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:)2 
2 

•TrT 


+  1      f    [     A     E(n(t1  )n(t?))    cos  2rfct1 
^2  ■'o  ^0  L 

cos   2irfct2  dti   dtp 

-  (±  it)2 

2 

+  L.  f   J        K^l't?)   aZ   cos  ^rrfpt-  (2.10) 

T2      0      0 

cos   2Tifct2  dt-    dt£. 

Substituting  equation  2,6  into  equation  2,10,  using  the 
properties  of  the  Dirac  delta  function  and  trignometric 
identities,  we  obtain 

E(rp)  =  (+  A_)   +  A2Nn  .         (2.11) 

*       2     Tr 
Therefore, 

Var(r?)  =  A2N0.  (2.12) 

Let  us  now  define  the  following  hypotheses: 

HQ:   r(t)  =  s0(t)  +  n(t),        (2.13a) 

Hji   r(t)  =  s^t)  +  n(t),        (2.13b) 

which  can  be  reduced  to 

H0:   r?  is  N(a£,  A2N0 ) ,         (2.14a) 
2    4T 

H-j   rP  is  N(-A2,  A2Nn ) ,        (2.14b) 
1  2   "Tfir 

where  N(m,n)  denotes  a  gaussian  random  variable  with  mean 

m  and  variance  n.   In  the  case  where  either  hypothesis 

is  equally  likely,  implying  that  either  a  binary  one  or 
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a  binary  zero  has  been  sent  with  equal  probability, 
and  where  the  penalty  for  making  an  incorrect  decision 
is  independent  of  whichever  hypothesis  is  true,  the 
decision  rule  for  minimum  probability  of  error  is 
always:  (see  for  example  3>M 

Choose  H0  if  r2  >  0,  (2.15a) 

Choose  E1    if  r2  <  0.  (2.15b) 

Thus,  we  have  shown  that  Figure  2.2  is  a  model  of 
the  optimum  receiver  for  our  signal  set  when  the  noise 
is  gaussian.  However,  we  have  not  yet  explained  why  we 
considered  only  additive  noise.   The  concept  of 
modelling  the  noise  as  inflicting  a  random  attenuation 
upon  the  information  part  of  the  received  signal  does  not 
change  the  optimum  receiver  structure  for  the  signal  set 
we  have  chosen.  All  operations  remain  the  same,  and  the 
decision  threshold  remains  at  the  origin;  however,  the 
probability  of  making  an  error  will  increase.   Since 
the  receiver  structure  remains  the  same,  we  feel  that 
there  is  no  necessity  for  considering  attenuation. 

However,  the  random  phase  shift  is  not  easy  to 
disregard.   The  primary  cause  of  random  phase  shift  noise 
is  the  fact  that  there  is  no  way  for  our  receiver  model 
to  insure  that  its  reference  oscillator  will  remain 
phase  locked  with  the  transmitter's  oscillator.   When 
unknown  phase  errors  exist,  the  receiver  loses  its 


-25- 
ability  to  distinguish  between  s,>(t)  and  s^(t).   The 
receivers  we  have  shown  in  Figures  2.1  and  2,2  are 
known  as  coherent  phase  receivers  since  they  assume 
that  their  oscillators  are  perfectly  locked  to  the 
transmitter  oscillator.   In  reality,  receivers  are 
phase  incoherent,  and  they  are  considerably  more 
complicated  than  the  simple  coherent  receivers  we  have 
shown,  (see  for  example  3»i+*5>6)  We  feel,  however, 
since  the  object  of  this  research  is  to  perform  a 
first  test  upon  the  dual  channel  receiver  design,  which 
we  present  in  the  next  section,  that  we  should  only  test 
the  simple  coherent  phase  receivers.  Moreover,  we  feel 
that,  since  both  receivers  we  are  comparing  will  be 
coherent  phase  types,  the  results  should  be  somewhat 
typical  of  the  results  to  be  expected  in  the  incoherent 
phase  case.  Finally,  if  our  results  indicate  signif- 
icantly Improved  performance,  then  the  more  difficult, 
and  more  realistic,  incoherent  phase  test  would  be 
justified. 

2.3  Joint  Channel  Receiver 

In  the  last  section,  we  presented  an  optimum 
receiver  for  use  when  the  corrupting  noise  is  gausslan. 
However,  when  the  noise  is  non-gaussian,  we  have  no 
guarantee  that  an  optimum  gaussian  receiver  will  yield 
optimum  performance.   With  knowledge  of  the  specific 
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non-gaussian  nature  of  the  actual  noise  encountered  in 
low  frequency  radio  communication  systems,  it  should  be 
possible  to  design  a  receiver  which  performs  better  than 
one  which  was  designed  for  gausslan  noise. 

Receiver  designs  which  have  been  tried,  or  proposed, 

for  dealing  with  the  non-gaussia.u  high  energy  noise 

1  2 

bursts  generally  Involve  some  type  of  nonlinear  device.  * 

These  receivers  use  clippers,  hole  punchers,  and  other 
similar  nonlinear  devices  to  limit  the  dominating  effect 

of  the  high  energy  noise  bursts. 

1 

Feldman  proposed  another  approach  to  the  receiver 

design  problem.  He  noted  that,  above  a  certain  threshold, 
the  noise  waveform  on  the  signal  channel  exhibited  a 
statistical  dependence  upon  the  value  of  the  noise 
envelope  on  a  nearby,  disjoint  bandlimited  frequency 
channel.  His  design  exploits  this  characteristic  in  a 
nearly  optimum  manner  by  applying  a  multiplicative 
weighting  function,  based  upon  the  inverse  of  the 
envelope  measurement,  to  the  signal  channel.   As  it 
applies  to  this  work,  this  receiver  design,  Figure  5-1 
in  reference  (1),  is  shown  here  as  Figure  2.3»   It 
should  be  noted  that  this  receiver  is  basically  a 
modified  linear  correlation  receiver. 

The  parameter  v  is  critical  to  this  design,  and 
it  is  given  by  equation  4.32  of  reference  (1),  which  is 
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Figure  2.3  Dual  Channel  Receiver 
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repeated  here: 

where  the  quantities  on  the  right  hand  side  are  given  in 
Table  1.1  of  this  work.   Whenever  the  noise  envelope  on 
the  second  channel,  v(t),  exceeds  the  threshold  value, 
vt,  there  exists  a  statistical  dependence  between  v(t) 
and  the  bandlimlted  noise  waveform,  y(t),  on  the 
information  channel.   Thus,  a  weighting  function,  z(t), 
given  below,  is  applied  to  the  information  channel: 
z(t)  =/~l,        v(t)^  vt 

[vt  /  v(t),  v(t)  >  vt.        (2.17) 
We  make  no  attempt  to  analytically  verify  that  this  in- 
verse first  order  weighting  function  is  indeed  nearly 
optimum.   Instead,  we  rely  upon  our  results  to  verify 
this  fact.   The  weighting  function  given  in  equation  2.17 
above  was  first  chosen  as  intuitively  pleasing. 

We  verified  our  choice  of  weighting  function 
experimentally.   First,  we  point  out  that  the  object  of 
the  dual  channel  receiver  is  to  reduce  the  problem  from 
one  in  which  we  have  a  constant  signal  power  level  and 
a  time  varying,  and  unknown,  (non-gaussian)  noise  power 
level  to  a  problem  in  which  the  noise  power  level  is  kept 
constant  and  the  signal  power  level  is  time  varying  in 
a  known  manner.  Keeping  this  objective  in  mind,  we 
conducted  simulations  for  all  four  weather/noise 
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condltions,  using  various  weighting  functions.   We 

simultaneously  simulated  an  ad  hoc  receiver  which  was 

composed  of  a  time  varying  signal,  formed  by  applying 

the  weighting  function  to  only  the  information  part  of 

the  received  signal,  and  a  gaussian  noise  process.   We 

then  chose  the  weighting  function  which  yielded  dual 

channel  receiver  performance  closest  to  the  performance 

of  the  ad  hoc  receiver.   The  final  weighting  function, 

z(t),  thus  determined  is  given  by  equation  2,17  above, 

2,4  Chapter  Summary 

In  this  chapter  we  have  presented  the  modem  to  be 
used  in  the  receiver  simulations  of  this  work.   In 
addition,  we  have  presented  an  analysis  of  a  linear 
correlation  receiver  which  operates  on  binary  antipodal 
signals  that  have  been  corrupted  by  additive  white 
gaussian  noise.   We  have  also  presented  an  improved 
design  for  a  low  frequency  digital  receiver,  which 
utilizes  information  about  the  specific  non-gaussian 
nature  of  the  low  frequency  atmospheric  noise  to  yield 
nearly  optimum  performance. 

In  the  next  chapter,  we  will  explain  some  of  the 
problems  encountered  in  attempting  a  digital  computer 
simulation  of  these  receivers.   We  will  also  give  examples 
of  the  calculations  necessary  to  relate  the  digital 
simulation  to  the  analog  system  being  evaluated. 


-30- 
Chapter  3 
SIMULATION  CALCULATIONS 

In  this  chapter,  we  will  present  the  information 
necessary  to  relate  our  digital  models  of  the  noise 
and  receivers  to  the  actual  analog  systems  we  simulated. 
In  the  first  section,  we  present  the  calculations 
necessary  to  convert  the  passband  receiver  model  to 
baseband.   In  the  second  section,  we  calculate  the 
probability  of  error  for  the  simulated  gaussian  noise 
case.   We  also  explain  how  we  checked  our  receiver  and 
noise  model  implementations  on  the  computer.   In 
Section  3»3$  we  give  the  calculations  which  relate  the 
sampled  data  noise  power  spectrum  to  the  analog  noise 
power  spectrum,  in  the  white  gaussian  noise  case. 

3.1  Baseband  Calculations 

The  received  signals  which  we  described  in  the 
last  chapter  were  located,  in  frequency,  in  the  receiver's 
passband.   It  is  a  well  known  fact  that  a  signal  must 
be  sampled  at  least  twice  as  fast  as  the  most  rapidly 
varying  frequency  component  in  the  signal,  in  order  to 
be  able  to  recover  the  signal  from  the  samples. 

The  highest  frequency  component  of  our  received 
signal,  r(t),  given  by  equation  3*1,    is  the  carrier 
frequency  plus  the  bandwidth.   The  received  signal,  with 
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y(t)  representing  the  bandlimited  atmospheric  noise,  is 
r(t)  =  +  A  COS  2rrfct  +  y(t).  (3.D 

In  our  case,  the  carrier  frequency  is  65  kHz,  and  the 
bandwidth  is  1  kHz;  thus,  the  highest  frequency  component 
allowed  in  r(t)  is  66  kHz,   Therefore,  to  accurately 
recover  r(t)  from  the  samples,  we  must  sample  the  signal 
at  a  sampling  frequency  of  at  least  132  kHz, 

'  Bit  lengths  for  low  frequency  digital  communication 
receivers,  the  length  of  time  in  which  the  receiver  must 
decide  which  of  the  two  possible  transmitted  signals  has 
been  received,  are  usually  about  20-50  msec.   For 
simulation  purposes,  we  have  chosen  the  lower  end  of 
this  scale,  in  order  to  keep  each  simulation  run  as 
short  as  possible.   Sampling  at  132  kHz  yields  a  time 
between  samples  of  about  0,0076  msec,  which  we  will  round 
up  to  0o01  msec  for  convenience.   To  simulate  one  bit  of 
information,  or  20  msec,  at  this  sampling  rate,  we  would 
require  2000  samples. 

Our  evaluation  of  receiver  performance  is  based 
strictly  upon  an  estimate  of  the  probability  of  an 
error  being  made  by  the  receiver  in  any  particular  bit 
of  binary  information  received.   In  order  to  form  this 
estimate,  we  have  chosen  the  following  approach:   count 
the  number  of  bits  simulated,  count  the  number  of  incorrect 
decisions  made  by  the  receiver,  and  divide  the  latter  by 
the  former,   01 sen  recommends  that  if  the  maximum  number 
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of  bits  we  simulate  is  1000,  then  the  minimum  probability 

of  error  we  can  estimate,  with  any  accuracy,  is  0.01. 
Further  details  about,  and  statistical  bounds  upon,  our 
estimation  procedure  will  be  given  in  Chapter  4. 

Based  on  the  above,  we  feel  that  1000  bits  is  the 
minimum  number  of  bits  we  must  simulate  in  order  to 
derive  meaningful  results.   In  this  case,  to  simulate 
the  passband  system,  we  will  need  2,000,000  samples  of 
the  analog  waveform,  taken  at  a  sampling  interval  of 
0.01  msec.   This  is  incredible,  especially  when  it  is 
realized  that  1000  bits  represent  only  20  seconds  of 
signalling.   This  is  also  impractical  for  various 
reasons:   the  first  reason  is  the  tremendous  computation 
time  needed,  and  the  second  is  that  our  noise  model  is 
incapable  of  accurately  simulating  such  fine  time 
structure  of  the  noise  process  as  is  required  by  this 
sampling  rate. 

We  have  thus  far  ignored  the  one  fact  which  allows 
us  to  solve  this  dilemma,  the  noise  and  transmitted 
signal  are  bandlimited.   Therefore,  we  can  accurately 
convert  the  passband  system  to  a  baseband  system  in 
which  the  maximum  frequency  component  is  1  kHz.   Clearly 
then,  the  maximum  sampling  rate  necessary  is  2  kHz,  and 
the  time  between  samples  is  0.5  msec.   Therefore,  the 
number  of  samples  per  bit  (20  msec  bit  length)  is  only 
40,  and  the  number  of  samples  per  1000  bit  simulation 
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run  Is  only  ^0,000.   Simulation  of  the  baseband  receiver 

is,  of  course,  a  much  more  tractable  problem. 

Any  bandlimited  process,  random  or  deterministic, 

can  be  reduced  to  lowpass  (baseband)  processes.   The 

procedure  for  doing  this  is  rather  simple,  and  it  is 

presented  in  various  references  on  communication  theory 

(see,  for  example  3»z+)»   We  simply  multiply  the  passband 

process  by  a  cosine  and/or  sine  of  the  carrier  (bandcenter) 

frequency,  and  pass  the  resulting  process  through  a  low 

pass  filter,   A  criterion  for  this  procedure  to  work  is 

that  the  bandwidth  of  the  passband  process  must  be  much 

less  than  the  carrier  frequency. 

Since  the  procedure  which  converts  a  passband  process 

to  a  baseband  process  is  actually  a  linear  filtering 

operation,  and  since  the  received  signal,  r(t),  is  a  linear 

combination  of  the  transmitted  signal,  s(t),  and  the 

corruptive  atmospheric  noise,  y(t),  we  can  treat  the 

conversion  of  each  of  these  processes  separately  and 

linearly  combine  the  results. 

First  we  consider  the  transmitted  signal  (we  consider 

only  s(t)  =  SQ(t),  since  s^(t)  follows  directly): 

s(t)  =  A  cos  27rfct.  (3.2) 

Then  the  baseband  process,  by  cosine  conversion,  s  (t), 

c 

before  low  pass  filtering,  is 

s    (t)   =  A  cos2(2TTfct)   =  A   (1   +   cos  47ifrt),      (3.3) 


and,  of  course,  after  filtering, 

sc(t)  -  A/2.  (3.^) 

Similarly,  by  sine  conversion,  we  get  s_(t): 

sg(t)   =  A   (cos   2TTfct)(sin   2trfct) 

=  A  sin  i+n-f   t,  (3.5) 

2  C 

and  after  low  pass  filtering,  thj.s  becomes 

ss(t)  =  0.  (3.6) 

Therefore,  we  see  that  in  order  to  convert  the  transmitted 

signal  portion  of  r(t)  to  a  baseband  process,  we  simply 

must  calculate  s_(t). 

c 

We  now  consider  the  baseband  conversion  of  y(t),  the 
bandlimited  low  frequency  atmospheric  noise  process.   By 
cosine  conversion,  before  low  pass  filtering,  we  obtain 
yc(t): 

yc(t)  =  y(t)  cos  27rfct.         (3.7) 
In  a  similar  manner,  before  low  pass  filtering,  ys(t)  is 

ys(t)  =  y(t)  sin  2irfct.  (3.8) 

After  low  pass  filtering,  y„(t)  would  be  known  as  the 
in  phase  component  of  the  baseband  noise,  and  y0(t) 
would  be  known  as  the  quadrature  component. 

At  this  point,  we  note  that  we  have  two  baseband 
processes  which  have  been  generated  from  the  bandpass 
received  signal  process,  r(t).   These  two  processes  are 
rc(t)  =  sc(t)  +  yc(t)  =  A/2  +  yc(t)     (3.9) 
and 

r.(t)  =  Sc;(t)  +  yAt)   =  yQ(t).       (3.10) 
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Since  r_(t)  contains  no  information  about  the  transmitted 

signal,  we  would  like  to  disregard  it  and  perform  a 
standard  baseband  correlation  receiver  operation  upon 
rc(t)  only.   It  can  be  shown-^  that  if  the  spectrum  of 
y(t)  is  even  over  the  bandwidth,  then  yc(t^)   and  ys(^2^ 
are  uncorrelated  for  all  observation  times  t^  and  t£» 
Further,  if  y(t)  were  a  gaussian  process,  then  yc(t^) 
and  ys(t2)  would  be  independent.   Then  rs(t)  would 
contain  no  additional  information  about  either  the  signal 
or  the  noise,  and  it  could  be  disregarded  with  no  loss 
in  optimality. 

Bandlimited  low  frequency  atmospheric  noise  is 

i 
definitely  non-gaussian,  but  Peldman  has  noted  that  it 

has  a  spectrum  which  is  determined  solely  by  the  band- 
limiting  filter's  frequency  response  squared,  which  we 
assume  to  be  even  over  the  bandwidth.   Therefore,  we  can 
state  that  yc(t^)  and  ys(t2)  are  indeed  uncorrelated  for 
all  time  instants  t^  and  t£«   Moreover,  we  assume  that 
the  statistical  dependence  between  yc(t^)  and  ys(to)  *s 
sufficiently  small  so  that  we  can  ignore  rs(t).   We 
feel  that  this  is  not  a  limiting  assumption,  especially 
when  viewed  in  the  light  that  it  allows  us  to  simulate 
a  baseband  receiver  as  opposed  to  the  nearly  impossible 
task  of  simulating  the  passband  receiver.   We  also  invoke 
the  fact  that  we  are  comparing  the  results  of  two 


-36- 
receivers  under  the  same  set  of  conditions  to  justify  our 

assumption  that  ignoring  rg(t)  is  of  little  consequence. 

Now  that  we  have  transformed  our  received  signal  into 
a  baseband  signal,  we  will  proceed  with  the  simulation 
of  a  baseband  receiver.   The  received  signal  is  now 
represented  by  equation  3«9»   As  a  practical  matter,  we 
generated  yc(t)  on  the  computer  by  sampling  y(t)  at  twice 
its  bandwidth  and  multiplying  the  samples  by  the  value 
of  cos  2irfct  at  the  sample  times.   In  the  next  section, 
we  will  state  how  we  tested  our  receiver  and  noise  model 
implementations  on  the  computer  by  calculating  the 
probability  of  error  in  the  gaussian  noise  case, 

3,2  Receiver  Implementations 

In  the  last  section,  we  described  the  manner  in 
which  we  converted  an  untractable  passband  receiver 
simulation  problem  into  a  more  tractable  baseband 
simulation.   In  this  section,  we  will  present  the  digital 
model  for  our  receivers.   We  will  verify  our  single 
channel  receiver  by  testing  it  with  gaussian  noise. 
We  will  check  the  simulation  results  with  analytically 
calculated  results,  and  we  will  show  that  these  results 
verify  our  digital  computer  implementation  of  the 
single  channel  receiver  as  correct.   We  will  then  test 
our  model  for  the  generation  of  the  atmospheric  noise 
samples  by  comparing  the  results  in  the  nearly  gaussian 
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quiet  noise  case  with  the  gausslan  results.   These  results 

will  be  similar  if  our  implementation  of  the  noise 

generator  is  correct.   Since  the  noise  generator  is  the 

same  for  all  noise  conditions  (quiet,  quiet  night>  tropical, 

and    frontal),  with  the  different  conditions  simulated  by 

changing  parameters  only,  we  will  then  assume  our  noise 

generator  is  correct  for  all  noise  conditions. 

The  basic  single  channel  correlation  receiver 

model  is  shown  in  computer  flow  graph  form  in  Figure  3*1. 

We  note  that  we  have  used  a  rather  simple  integration 

algorithm,  x^hich  is  a  slightly  modified  trapezoidal 

rule.   We  explain  it  as  follows: 
rT  N 

1  J   x(t)  dt  ~1  g  (x(n)  +  x(n-l))  /t,   (3.11) 
TO  T  n=l  2 

where  it  is  the  Interval  between  samples,  NAt  =  T,  and 

x(n)  is  a  sample  of  x(t)  at  time  nAt.   We  note  that,  in 

our  case,  we  are  Integrating  over  one  bit  length,  and 

we  can  only  include  N  samples  in  our  summation.   This 

is  because  the  sample  called  x(N)  for  the  i   bit  is 

th 
also  the  sample  called  x(0)  for  the  (i+1)   bit.   Since 

the  signal  on  each  of  these  bits  may  be  different,  we 

must  assign  the  end  point  to  one  bit  only.   We  therefore 

consider  only  samples  x(0)  through  x(N-l)  for  each  bit, 

and,  in  order  to  simplify  our  integrator,  we  assume  that 

x(N)  would  approximately  equal  x(0).   Thus, 

N-  N-i 

1  S.  (x(n)  +  x(n-l))  At  ^  1  <T  x(n),    (3.12) 

T  n=l  2  ~  N  n=0 


-38- 


(,,  -•  o  -  tf   jficorret*- 

ir-o   iil    of  W.+ 


^__*r 


t^f  t     Si^al    -  A 


£  -  0 


ik £ 


-t  -     k  Ai      +  jv  t, 


Figure  3,1  Single  Channel  Receiver 
Computer  Plow  Graph 
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which  makes  our   integrator 

rT  N-l 

1    /       x(t)   dt-1    <i    x(n),  (3.13) 

TO  N  n=0 

which  is  the  slightly  modified  trapezoidal  integration 
rule  that  we  have  used.   We  feel  that  when  N  is  relatively- 
large  (for  example,  40,  as  in  our  case),  the  errors 
introduced  by  this  simplification  are  small. 

In  Figure  3»2,  we  show  the  results  of  a  simulation 
of  our  single  channel  in  which  the  noise  was  a  gaussian 
process.   We  conducted  three  simulation  runs  of  1000  bits, 
at  each  of  two  bit  lengths,  averaged  the  probability  of 
error  estimates  for  each  bit  length,  and  plotted  them 
versus  the  signal-to-noise  ratio  (SNR).   We  also 
calculated  the  probability  of  error  for  each  SNR  simulated, 
and  those  curves  are  plotted  along  with  the  corresponding 
simulation  results.   As  can  be  seen,  the  calculated  and 
experimental  curves  are  in  close  agreement.   We  therefore 
conclude  that  our  digital  computer  implementation  of  the 
single  channel  receiver  is  correct. 

In  Figure  4.4,  we  present  the  results  of  our 
simulation  runs  using  our  atmospheric  noise  generator  with 
the  quiet  noise  parameters.   We  note  that,  as  we  expected, 
the  results  in  this  case  are  nearly  the  same  as  in  the 
gaussian  noise  case,  for  the  single  channel  receiver. 
We  feel  that  these  results  verify  our  implementation  of 
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Figure  3.2  Gaussian  Results 
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the  atmospheric  noise  model  in  the  quiet  noise  case,  and 
hence  in  all  other  cases,  since  they  involve  simple 
parameter  changes  in  the  generator. 

Finally,  Figure  3»3  shows  a  computer  flow  graph  of 
the  dual  channel  receiver.   We  see  that  it  is  basically 
the  same  as  the  single  channel  implementation,  with  the 
exception  that  it  incorporates  the  weighting  function, 
which  is  based  upon  the  second  channel  noise  envelope 
samples,  inside  the  integration  routine.   We  feel  that, 
since  our  results,  shown  in  Figure  4.1  through  Figure 
4,4,  indicate  significant  improvements  in  performance  in 
the  more  non-gaussian  noise  cases  and  comparable  performance 
in  the  nearly  gaussian  quiet  noise  case,  when  compared 
with  the  single  channel  receiver,  our  dual  channel 
receiver  implementation  is  correct, 

3,3  Analog  to  Digital  Calculations 

In  testing  our  single  channel  receiver  versus 
analytically  calculable  results,  we  encountered  the 
problem  of  relating  the  analog  and  digital  noise  power 
spectra.   We  will  calculate  that  relationship  in  this 
section.   First,  however,  we  will  show  how  we  calculated 

the  probability  of  error  per  binary  bit  of  information 

3 
when  the  noise  is  gaussian. 
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Define  the  average  power  in  the  information  part  of 

the  received  signal,  r(t),  as  follows: 

Ps  =  Es/T,  (3.14) 

where  T  is  the  bit  length  and 


=  J       s*(t)  dt,  (3.15) 


Es 
S   ^  0 

where  s(t)  is  the  baseband  signal  portion  of  r(t),  i.e. 

s(t)  =  +  A/2.   Define  the  source  rate,  R,  as  the  rate 

at  which  binary  bits  are  generated  by  the  transmitter. 

Then,  the  energy  per  bit,  Eb,  is  given  by 

Eb  =  Ps/R. 

Therefore, 

E,  =  1/T  /   s2(t)  dt  =  (A/2)".     (3.16) 

D      y0 R 

R 

The  baseband  SNR  of  the  simulated  signal  and  noise 

is  defined  as  follows: 

SNR  =  Eb/NQ,  (3.17) 

where  NQ/2  is  the  variance  of  each  independent  sample  of 

the  gaussian  noise  process  in  the  digital  computer 

simulation.   However,  in  order  to  analytically  compute 

the  probability  of  error,  we  must  obtain  an  expression 

for  the  power  spectrum  of  a  continuous  white  gaussian 

noise  process,  which  is  equivalent  to  the  sampled  noise 

process  used  in  the  simulation.   We  can  do  this  in  the 

following  manner,  in  which  n(.)  represents  the  noise 


1 


process: 

fcVi 

.   the  k        noise   sampl 

■At 


Let  the  kth  noise  sample  be  defined  as 


n(k)  -  1    J         n(t)   dt,  (3.18) 

2Jt    o 


then 

-At  ^At 


E(n?(k))   =  N0/2  =   E(l J  n(tt)   n(t?)   dtj   dt2) 

=  J__  J      J  E(n(ti)   n(t2))   dtj   dt2.      (3.19) 

At2      ° 

Recognizing  that 


At2      0     "0 


E(n(t1)  n(t2))  t  Neq  Jd^-tg),       (3.20) 

where  £(.)  is  the  Dirac  delta  function,  we  get 

NQ/2  =  NeQ/2,  (3.21) 

At 
and  therefore, 

Neq  =  %Atf  (3.22) 

where  At  is  the  time  interval  between  samples  and  N  „  is 

r  eq 

the  magnitude  of  the  noise  power  spectrum  of  a  continuous 
white  gaussian  process,  analogous  to  the  simulated  process. 
Thus,  the  analog  SNR  corresponding  to  the  digital  SNR 
(equation  3,17)  is  given  by 

SNR  =  Eb/Neq  =  Et/NgAt.  (3-23) 

For  antipodal  signalling  schemes,  the  probability 

of  error  per  bit,  with  gaussian  noise,  can  be  calculated 

3 
as  follows: 

P(E)  =  Q(//~2Eb/Neq'),  (3.24) 

where  Q(x)  is  defined  as 


Q(x)  =  J   1_  e"u  /2  du,  (3.25) 

x  2tt 

whose  complement,  1  -  Q(x),  is  tabulated  in  various 

references  (e.g.  8).   Thus,  in  terms  of  the  parameters 

of  our  digital  simulation,  the  probability  of  error  per 

binary  bit,  with  white  gaussian  noise,  can  be  analytically 

calculated  by  means  of  the  Q( . )  function  as 

P(E)  =  Q(/V2Eb/N0At' ).  (3.26) 

3.^-  Chapter  Summary 

In  this  chapter,  we  have  explained  why  and  how  we 
converted  our  bandlimited  receiver  from  a  passband  to  a 
baseband  system  for  digital  simulation  purposes.   We  went 
on  to  explain  how-  we  checked  our  implementations  of  the 
receiver  and  noise  models  by  comparing  simulated  and 
analytical  results  using  gaussian  noise  corruption  and 
corruption  by  the  quiet  noise  case  of  atmospheric  noise. 
Finally,  we  gave  the  relationship  between  the  simulated 
noise  power  spectrum  and  the  analogous  continuous  power 
spectrum,  which  was  necessary  to  calculate  the  probability 
of  error  per  binary  bit  analytically. 

In  the  next  chapter,  we  will  present  the  results  of 
our  simulation  runs  using  the  atmospheric  noise  model 
described  in  Chapter  1.   We  will  show  the  averaged 
results  for  five  runs  of  each  of  the  four  basic 
weather/noise  conditions:   quiet,  quiet  night,  tropical, 
and  frontal.   We  will  show  how  we  obtained  our  estimate 
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of  P(E)  per  bit,  and  we  will  give  bounds  on  that  estimate. 
We  will  also  give  an  example  of  the  problems  we 
encountered  as  we  proceeded  with  the  simulations.   Finally, 
we  will  summarize  the  conditions  under  which  our  results 
were  obtained  so  that  the  reader  may  evaluate  their 
usefulness. 
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Chapter  4 

SIMULATION  RESULTS 

Before  presenting  our  results,  we  will  briefly 
summarize  our  work  up  to  this  point.   We  discussed,  in 
Chapter  1,  bandlimited  low  frequency  atmospheric  noise, 

its  causes,  its  non-gaussian  nature,  and  a  model  for 

i 
this  noise.   The  model,  developed  by  Feldman,   is 

especially  suited  for  simulation  on  a  digital  computer. 

In  Chapter  2,  we  discussed  several  receiver  schemes, 

which  could  be  used  to  decide  which  of  tx-ro  antipodal 

signals,  corrupted  by  atmospheric  noise,  was  transmitted. 

Specifically,  we  discussed  a  linear  correlation  receiver, 

which  is  optimum  when  the  noise  encountered  is  gaussian. 

We  then  discussed  a  design,  based  upon  the  observed 

properties  of  the  non-gaussian  noise,  which  is  considerably 

more  complicated  than  the  simple  correlation  receiver, 

but  which  should  yield  optimum  performance  in  the 

presence  of  this  noise.   Finally  in  Chapter  3j  we 

discussed  the  problems  involved  with  a  simulation  of 

these  receivers  and  bandlimited  low  frequency  atmospheric 

noise. 

In  this  chapter,  we  present  the  results  we  obtained 

in  our  simulation  of  these  receivers.  We  tested  each 

receiver  for  each  of  four  weather/noise  conditions, 

which  our  noise  model  is  capable  of  simulating,  ranging 
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from  the  nearly  gaussian  quiet  noise  case  to  the  extremely 
non-gaussian  frontal  noise  case.   For  each  receiver  and 
weather/noise  condition,  we  made  five  simulation  runs, 
each  run  simulating  the  reception  of  1000  bits  of  binary 
information  by  the  receiver.   In  order  to  characterize 
receiver  performance,  we  formed  an  estimate  of  the 
probability  of  error  per  bit  of  received  information  for 
each  run.  We  then  averaged  the  results  of  the  five  runs 
for  each  case  and  plotted  these  averages  as  functions  of 
the  signal-to-noise  ratio  of  the  received  baseband 
signals. 

We  present  a  discussion  of  the  method  in  which  we 
formed  our  estimates  of  probability  of  error,  and  we 
give  bounds  upon  our  estimating  procedure.   We  also 
bound  our  results  by  presenting  tables  of  certain 
statistics  of  the  results  of  the  five  individual  runs 
for  each  case.   We  discuss  briefly  one  problem  we 
encountered  during  the  simulation  runs,  and  why  this 
problem  existed  and  how  we  solved  it.   Finally,  in  the 
last  section,  we  present  a  summary  of  this  research, 

4.1  Results  and  Bounds 

In  Figures  4,1  through  4.4,  we  have  plotted  the 
results  of  our  simulation  runs  for  frontal,  tropical, 
quiet  night,  and  quiet  weather/noise  conditions, 
respectively.   The  curves  shown  are  plotted  on  semi-log 
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paper,  with  the  ordinate  representing  the  probability  of 

error  per  bit  and  the  abcissa  being  the  signal-to-noise 
ratio  in  decibels  (20  log  SNR).   There  are  three  curves 
shown  in  each  figure.   These  curves  represent  the  results 
of  the  single  channel  correlation  reciever,  the  dual 
channel  receiver,  and  the  ad  hoc  receiver  (described  in 
Section  2.3)  which  we  used  as  a  check  upon  our  dual 
channel  receiver. 

The  curves  shown  in  each  graph  represent  the  average 
of  five  simulation  runs  of  1000  bits  of  received  infor- 
mation each.   The  bit  length  for  all  runs  was  20  msec, 
making  the  simulated  length  of  each  run  20  seconds.   The 
signal-to-noise  ratio  scale  was  defined  in  the  following 
manner:   It  is  necessary  to  realize  that  the  noise  power 
level  is  vastly  different  from  one  weather/noise  condition 
to  any  other.   However,  the  banc? limited  low  frequency 
atmospheric  noise  model  used  in  these  simulations  consists 
of  two  linearly  combined  functions,  a  background  low  level 

gaussian  process,  and  a  process  which  represents  the  high 

1 
energy  burst  behavior  of  the  noise.    The  background 

gaussian  noise  is  assumed  to  have  been  generated  in  the 

2 
receiver  itself  (receiver  front-end  noise),   and  therefore, 

it  is  independent  of  the  particular  weather/noise  condition 

being  simulated.   Since  we  would  like  our  SNR  scale  to 

be  independent  of  the  weather/noise  condition  being 
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simulated,  we  base  it  upon  the  signal  power  and  the  power 
level  of  the  background  gaussian  process. 

Thus,  our  SNR  is  given  by  equation  3.23,  which  is 
repeated  here: 

SNR  =  Eb/N04t, 
where  Eb  is  the  energy  per  bit,  as  given  by  equations  ^*lk 
through  3«l6,  ^q/2   is  the  variance  of  each  independent 
sample  of  the  background  gaussian  noise  process  of  our 


noise  model  (all  the  parameters  of  our  noise  model  have 

!0< 


1 
been  normalized  such  that  Nn/2  ~  1),   and  A t  is  the 


sampling  interval. 

The  estimate  of  the  probability  of  error  per  bit- 
was  formed  by  counting  the  number  of  incorrect  decisions 
made  by  the  receiver  and  dividing  that  number  by  the 
total  number  of  bits  simulated.  For  example,  in  a  run 
where  we  simulated  1000  bits,  and  where  the  receiver 
made  351  incorrect  decisions  (we  have  assumed  that  the 
probability  of  making  an  error  is  independent  of  the 
signal,  either  a  binary  one  or  a  binary  zero,  which  was 
sent),  our  estimate  of  the  probability  of  error  per  bit 
in  that  run  would  be  0,351, 

We  now  present  an  evaluation  of  our  estimation 
procedure.   First,  we  assume  that  the  probability  that 
an  error  is  made  on  any  particular  bit,  for  any 
particular  receiver  and  weather/noise  condition,  is 
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constant  and  can  be  denoted  by  the  value  p.   Then  the 
probability  that  a  correct  decision  is  made  on  any  bit 
is  simply  (1  -  p)#   We  assume  that  the  decision  made  on 
each  bit  is  independent  of  any,  and  all,  other  bits. 

Now  we  can  assign  a  random  variable,  say  k*  ,  to  the 
outcome  of  the  i   decision  of  the  receiver.   The  random 
variable,  k  ,  is  defined  to  have  the  following  probability 
density  function? 

Pd^  =  0)  =  1  -  p,  (4.1a) 

P(ki  =  1)  =  p.  (4.1b) 

Therefore,  we  see  that  k.  is  one  if,  and  only  if,  an 
error  is  made  on  the  i   bit,  and  k*  is  zero  otherwise. 
Now  we  define  the  random  variable  s,  such  that 

N 
s  =  1/N  ^  k,,  (4.2) 

i=l 

where  N  is  the  number  of  bits  simulated.   We  know  that 

the  set  of  k,  is  a  set  of  independent  identically 

distributed  random  variables,  each  having  the  probability 

density  function  given  in  equations  4,1  above,  which  we 

note  is  the  well  defined  binomial  density.   Therefore, 

o 

the  mean  and  variance  of  each  k^  are  given  below: 

E(ki)  =  p,  (4.3) 

Var(k1)  =  p  (1  -  p) 

=  P  -  P2.  (4.4) 

Now,  we  see  that,  by  the  central  limit  theorem,  if  N  is 
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large,  the  probability  density  function  for  s  is  gaussian 

with  mean  and  variance  described  by 

N 
E(s)  =  E(l/N  ^  k,  ) 
1=1   x 

N 

=  1/N  ^L   E<ki  ) 
1=1     x 

=  p,  (4.5) 

and 

Var(s)  =  E(s2)  -  (E(s))2 

N   N 
=  E(l/N2  ;>   ^  k,  kj  -  p2 
1=1  j=l     J 

-  (P  -  P2)/N,  (4.6) 

where  we  have  used  the  independence  of  the  Vl**      We  make 

note  of  the  facts'  that  s  is  an  unbiased  estimate  of  p, 

and  that  the  variance  of  s  decreases  toward  zero  as  the 

number  of  bits  simulated  increases,. 

Finally,  in  order  to  establish  a  bound  upon  our 

estimates  of  probability  of  error  per  bit,  as  given  by 

p 
s  in  equation  4.2,  we  used  a  table  of  the  complement 

of  the  Q(.)  function,  described  in  equation  3«25»   We 

note  from  the  table  that,  with  95%   confidence,  our 

estimate  of  the  probability  of  error  per  bit  will  be 

within  two  standard  deviations  of  the  actual  value.   We 

also  note  that  the  variance  of  our  estimate,  and  hence 

its  standard  deviation,  is  greatest  when  the  actual 

value  of  p  is  0.5*   We  therefore  use  this  case  for  a 
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worst  case  analysis  of  our  estimation  procedure. 

When  p  =  0.5»  the  variance  of  s  is 

Var(s)  =  (0,5  -  0.25)/N  -  0.25/N.      (4.7) 
In  our  case,  we  have  made  five  runs  of  1000  bits  each; 
thus  making  N  =  5000.   Therefore, 

Var(s)   =   0.25/5000,  (4.8) 

and  the   standard  deviation,  <Tst    is 

(Ts  =  J  0.25/5000'  =  0.00708,  (4.9) 

and   thus, 

2<TS  =   0.01416.  (4.10) 

From  the  above,  we  can  see  that,  with  95%   confidence, 
our  estimates  of  the  actual  probability  of  error  per 
binary  bit  should  be  within  less  than  +  0.014  of  the 
actual  values.   As  a  matter  of  practical  concern,  we  note 
that  each  run  of  1000  simulated  bits  of  information  was 
conducted  independently  of  the  others.   We  therefore 
also  present  a  summary  of  the  results  of  the  five  runs 
for  each  weather/noise  condition  and  receiver.   These 
results  are  presented  in  tabular  form  in  Tables  4.1 
through  4.4.   In  these  tables,  we  give  the  average  of 
probabilities  of  error  of  the  five  runs  for  each  case  and 
the  maximum  deviation  from  that  average. 

4.2  Difficulty 

The  most  serious  difficulty  encountered  in  this 
research  was  an  apparent  breakdown  in  our  single  channel 


-58- 


Sio^le   cn«-«ne\ 

c^uko.l    cWimtl 

Oc{       'lUC- 

r 

t'cO 

-Z8 

nc   A^fr- 

nc    d  a^- ci- 

Ml- i  a  i  a. 

-%X 

uc    cS-Ct-to- 

ne    <&a"t  a. 

He   Ja(<*- 

~\y 

He     rioftc*- 

vie    c^a-tcL 

no    Jia-t  *- 

-(0 

ho    &ut<^ 

lac     Jus-t  £<- 

no    c^a  t  &* 

-n 

O.U,$$               0,022. 

C.3li*           c,c;5 

c,n°\        ceo*) 

*  * 

9 

O  ■  MSlr                 CO  13 

C.37C            CCC'i 

c,334          0.023 

O.  MS3              0.<?M 

c33o         cccci 

C  2-8  1              0.0/5 

(4 

©.    U  "7 "?                CC/I 

C/XIS            O.OlZ 

0.2  15              C.f/i 

2c 

0      U(rc'\                COM 

c.234            O.OlS 

C./45"           cci] 

56? 

c.  H5l          o■o^'b 

an          cci'i 

con          cool 

.3* 

0  .  450              O.O20, 

c.  1 1  9           c  -  oi  ? 

C.CX9,            0.CC5 

3$ 

c  .  43^            cc^3 

C.ClS             O.Ciq 

o.ooi         o.coS 

^ 

© ,  «4i a             C.  CX) 

ccm         0.007 

e.ooi         o.ooi 

50 

o,3^           o.  wi 

C-C2X               C.CCZ 

O  ,  CCO             COCO 

Table  h%\   Frontal  Noise 


.<>  ivM^le  cV, 

Xftflt  \ 

dluco  Gnavv 

ne\ 

Ad      Wco 

.SMn(^)    ' 

Wei 

i>m,i(.  df\T.   j  Jul  ra.c,(i 

mat-  a^,   a  uera  h  t 

mat,  dfiT, 

-2$ 

mm — 

(T.USi 

C0.CIC           CMll 

0,C2.f    i     O.^l1) 

0,0  a* 

-71 

<T.47t 

C   -C\C 

0,3ci3 

O.023       £,3^1 

c.c^i 

-J<a 

C  Mfc5 

O.on 

C3S3 

0,0/7 

c?,3M5 

CC2Z 

-  \C 

c.  aq^) 

0  >0%\ 

0-^1 

0<  025  |    0,2?^ 

c5.r/3 

-  4 

C-  tUS 

0.03S 

C.22C 

C  ,011  |    0 ,  2. 1  b 

0,C2^ 

+-  1 

6.   L|ol 

C.D^S 

C  .  I  4 c/ 

o.oai  !  031 

0,017 

% 

0.374 

C.Cqx 

0,077 

0.  C  »8     0.062 

0.013 

14 

0.  341 

c.053 

0.03k 

CCI/    U> .  0  /  Cs 

0.01  X 

ac 

0 ,  3  C  (* 

C ,  c  5G 

0.  0/7 

<5.(TC3    |  0,  003 

C  Do3 

at. 

O,  £70 

0  .c$2 

o.ooi 

C7.00  4        0.  00 0 

o.coo 

33 

0.2.3) 

d.c>44 

O-OOL) 

O.cc  3    !  c.  ccc 

c.ccc 

33 

o.  j?3 

0.034 

0.00% 

C  .CO  7      1   0.  CCO 

C.cc  0 

44 

© .  1 1\  0 

o.  'i1  5l 

0-00) 

0 , : ;  /     j  c ,  c  c  c 

o.ooo 

£0 

C,c%l 

0  ■  c  /3 

o.coo 

C.ccc      C.cco 

o.ooo 

» 

\ 

Table  km2   Tropical  Noise 


5vvk.\c    (.V 

0  vir  I  f"  V 

c^ufA  o\ 

i  unc  1 

cA  L- < 

-59  - 

sueoi 

P    pt6> 

n\a%,  (it  >t 

(1  fe> 

,   if-v 

iuft*.  Ol?l! 

.    ~ag 

0.4fcZ 

G.C3! 

0,407 

ecu 

0.i|W 

0.033 

.       ~*2 

O.MMG, 

c.cafc 

0.37g 

Cl'H 

>TJ 

C.035 

-IU 

o.4a3 

6.C3S 

0-332 

0.  022 

)  4S 

cell 

-  10 

C.3S3 

0.03? 

<?.^gc 

0.C  (7 

C.l°\  % 

c  .cxi 

-M 

0.3^^ 

6-C3C 

0.2/4 

0.019 

c.i\$ 

C  .0X1 

t  2 

0.3C^ 

0.039 

0,(37 

ecu 

0.  )4fe 

O  ,0  2.0 

8 

o.*m 

0.040 

C?.067 

0.01b 

6' .  6  7  / 

0  .cxo 

H 

0,  J  83 

C.C3S 

0,62? 

0  .on 

o.ou 

C,c\  | 

2.0 

0,f3M 

0.031 

C.COS 

O  .010 

0 .  o  o  3 

Coot 

2G 

0,0$? 

O.CiH 

CCCX 

0.003 

coco 

O.ccc 

3* 

0 . 0  5  2 

CCll 

Coco 

C.CCI 

C .  0  CQ 

o.  000 

3* 

0.076 

©,Co8 

c.ccc 

c.ccc 

c.ccc 

O.GOO 

MM 

0.oQ$ 

o.cot) 

o.ccc 

CcCO 

0,  COO 

0,  CCC 

£"0 

CCC  1 

ccoi 

0  .  000 

o.o  co 

o.ccc 

c.  ooo 

Table  4 

.3  Quie^ 

t  Night 

Noise 

txck  hoc 

■ 

rrik-t.&e-C. 

oMi7c.de< 

~28 

0.333 

O,0«lfe|  6.3c;r 

0  iCZC 

) 

O.  MC'4_ 

0.07C 

-22 

0.34<© 

coze 

c.3>3 

C.cz  | 

©•374 

0,  0^ 

-  tfe 

o.2«:s 

O.023 

caW 

0.023 

0332 

0.022 

-\c 

0.21b 

0  •  02fc, 

o.?32 

o.oa'fc 

6.2.7C 

o  ■  o  x  c 

-  U 

0.I3& 

0,0X1 

0'.  lb  5" 

c.c27 

0.  11*1 

0     OXQ, 

+  2 

0.05? 

com 

o-ciS 

CCll 

0>- 114 

o.otl 

? 

0.0(2 

C-.ocM 

o.oi? 

O.CCG 

0-0)1 

con 

m 

0.001 

CD.ooi 

0.002 

o,ooi 

O.OOio 

0.OCL 

.    ?-° 

O  ■  ceo 

COCO 

o .  coo 

C.CCO 

o.ooo 

ccoi 

aw 

C.  OC  0 

£>.OO0 

0  .  c  0  0 

C.CCO 

O.ccc 

o.ccc 

32 

fitOOO 

0.000 

o.ooc 

G .  C  00 

o.ooo 

o.ooo 

3* 

ex  ooo 

o.ooc 

C  .CCC 

O.OCG 

0.  000 

o.oo  o 

44 

a  ooo 

O.  OC  0 

0,660 

O.O0C 

o  ooo 

O.OOO 

£0 

a  ooo 

0,000 

o.ccc 

O-OOO 

o.coc 

0.  OC  0 

Table  4.4  Quiet  Noise 


-60  - 
linear  correlation  receiver  model.   Our  original  intention 
was  to  make  only  one  simulation  run,  of  1000  bits,  for 
each  receiver  and  weather/noise  condition  while  varying 
both  signal  magnitude  and  bit  length.   We  note  that  an 
increase  in  either  signal  magnitude  or  bit  length  results 
in  an  increase  in  SNR,  and  therefore  should  result  in  a 
decrease  in  the  probability  of  error  per  bit.   For  our 
dual  channel  receiver,  this  was  indeed  the  case,  but  with 
the  single  channel  model,  we  found  that  the  probability 
of  error  was  increasing  slightly  with  increasing  bit 
length  for  all  signal  magnitudes  below  a  certain,  non- 
constant  threshold. 

We  reasoned  that  it  would  be  interesting  to  see  how 
well  each  receiver  performed  as  a  function  of  both  the 
signal  magnitude  and  the  bit  length.   Therefore,  we 
thought  that  we  were  keeping  all  other  factors  constant 
and  varying  only  the  factor  of  interest,  either  the  bit 
length  or  the  signal  magnitude,  from  run  to  run.   We 
must  diverge  here  slightly  and  explain  the  manner  in 
which  we  generate  the  noise  process  upon  a  digital 
computer.   To  generate  any  random  process  upon  a  digital 
computer,  we  in  fact  are  generating  a  sequence  of 
pseudorandom  numbers,  these  numbers  being  samples  of  the 
noise  process.   In  order  to  generate  this  sequence,  we 
must  initialize  all  random  number  generators  in  the 
computer  program  (see  Figure  1,1).   Furthermore,  if  we 
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initialize  all  the  random  number  generators  the  same  way 

in  every  rim  of  the  same  program,  we  will  generate  the 

o 
identical  random  number  sequence. 

Now  we  can  see  that  when  we  varied  the  signal 
magnitude  only,  and  kept  the  bit  length  constant,  the 
noise  process  which  corrupted  each  bit  was  the  same  for 
all  signal  magnitudes  since  we  initialized  all  random 
numbers  in  the  same  manner  at  the  start  of  each  run. 
However,  when  we  varied  the  bit  length,  while  keeping 
the  signal  magnitude  constant,  the  noise  process  seen  by 
each  bit  changed  radically  with  bit  length,  regardless 
of  how  we  initialized  the  random  number  generators.   Thus, 
we  were  varying  two  factors,  both  the  noise  process  and 
the  bit  length.   We  made  several  more  runs,  with  different 
initial  conditions  upon  the  random  number  generators,  and 
we  found  that  the  probability  of  error  did  indeed  decrease 
as  we  increased  bit  length  in  these  runs.   We  therefore 
reasoned  that  our  first  run  had  been  made  with  a  non- 
typical  noise  process  that  caused  the  probability  of 
error  to  increase  slightly  with  Increasing  bit  length. 

We  subsequently  decided  that  we  should  make  several 
runs,  starting  with  different  initial  conditions  upon 
the  random  number  generators  each  time,  for  each  bit 
length  and  each  signal  magnitude0   However,  due  to  time 
constraints,  we  found  that  we  could  only  make  several 
runs  for  each  weather/noise  condition,  receiver,  and 
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signal  magnitude  at  one  bit  length.  We  feel  that  in  this 
way  our  results  are  typical  of  the  actual  results  to  be 
expected  for  that  bit  length.  We  chose  a  bit  length  of 
20  msec,  about  the  shortest  bit  length  practical  in  low 
frequency  digital  communication  systems,  because  it 
resulted  in  the  shortest  time  per  simulation  run, 

4,3  Summary  of  Research  Conditions 

In  this  section,  we  present  a  summary  of  the 
receiver  and  noise  models  we  used  in  the  simulation 
runs  which  we  made.   We  also  give  a  short  summary  of 
the  results  we  obtained. 

The  model  for  bandlimited  low  frequency  atmospheric 

1 
noise  which  we  used  was  developed  by  Feldman,    This  noise 

model  accurately  models  the  first  order  statistics  and 

the  time  structure  of  the  actual  bandlimited  noise 

waveform  for  four  weather/noise  conditions.   These  four 

conditions  include  the  nearly  gaussian  quiet  noise  case, 

the  quiet  night  noise  case,  the  tropical  noise  case,  and 

the  most  severe  frontal  noise  case.   These  four  conditions 

model  noise  which  is  typical  of  large  of  areas  of  the 

world  over  significantly  large  periods  of  time.   The 

quiet  noise  conditions  exist  when  the  atmsophere  in  the 

vicinity  of  the  receiver  consists  of  a  stable  air  mass. 

The  tropical  conditions  exist  when  the  air  is  warm,  moist, 

and  unstable.   Quiet  night  conditions  represent  a 
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transitional  period  between  quiet  and  tropical  conditions. 
And,  of  course,  frontal  noise  conditions,  the  most  violent 
conditions,  exist  when  there  is  a  front,  with  much 
thunderstorm  activity,  in  the  vicinity  of  the  receiver. 
A  canonic  noise  generator,  useful  in  computer  simulations, 
and  a  table  of  associated  parameters  are  given  in  Figure 
1.1  and  Table  1.1  of  this  work. 

For  our  simulation,  we  chose  binary  antipodal 
signalling,  with  a  coherent  PSK  modem.   The  carrier 
frequency  of  our  signals  was  6$   kHz  and  the  receiver's 
bandlimiting  filters  had  a  bandwidth  of  1  kHz.   We  also 
measured  the  noise  envelope  on  a  second  frequency  channel, 
with  bandwidth  1  kHz,  centered  at  83  kHz.   These  frequencies 
were  chosen  because  they  are  typical  of  digital  communication 

systems,  and  we  have  the  parameters  for  simulation  of 

1 
the  noise  model  at  these  frequencies. 

As  a  standard  of  comparison,  we  simulated  a  single 

channel  linear  correlation  receiver,  using  a  20  msec  bit 

length,  operating  in  the  above  described  low  frequency 

atmospheric  noise.   Our  test  statistic,  an  estimate  of 

the  receiver's  probability  of  error  per  binary  bit,  was 

calculated  simply  by  counting  the  number  of  incorrect 

decisions  made  by  the  receiver  and  dividing  by  the 

number  of  bits  simulated.   We  made  five  runs,  of  1000 

bits  each,  for  each  weather/noise  condition  for  various 

slgnal-to-noise  ratios. 
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We  then  tested,  in  a  similar  manner,  a  dual  channel 
receiver  which  uses  the  fact  that,  above  a  certain 
threshold,  the  value  of  the  noise  envelope  on  the  83  kHz 
channel  and  the  value  of  the  noise  waveform  on  the  signal 
channel  are  statistically  dependent.   Due  to  the  sampling 
rates  involved,  and  limitations  of  the  noise  model,  we 
simulated  both  of  the  above  receivers  in  the  baseband 
rather  than  the  passband. 

We  found  that  this  dual  channel  receiver  performed 
better  than  the  single  channel  receiver  for  the  three 
worst  noise  conditions  (frontal,  tropical  and  quiet  night), 
and  as  the  signal-to-noise  ratio  increased,  the  performance 
improvement  exceeded  an  order  of  magnitude.   However,  in 
the  quiet  noise  case,  we  found  that  the  dual  channel 
receiver  yielded  approximately  the  same  performance  as 
the  single  channel  model, 

^4-, 4  Chapter  Summary 

In  this  chapter,  we  have  presented  our  results  of 
the  simulations  of  the  single  and  dual  channel  receivers. 
We  also  presented  an  analysis  of  the  estimate  of  probability 
of  error  per  bit  which  we  used,  and  we  gave  bounds  upon 
our  results.   We  explained  a  difficulty  we  encountered 
in  the  course  of  the  simulations,  and  we  explained  how 
we  overcame  it.   Finally,  we  presented  a  summary  of  the 
research  reported  in  this  work.   In  the  next  chapter,  we 
present  our  conclusions  and  ideas  for  further  research. 
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Chapter  5 
CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

We  first  present  a  brief  summary  of  this  research. 
In  Chapter  1,  we  described  low  frequency  atmospheric 
noise  and  a  model  for  this  noise.   In  Chapter  2,  we 
presented  two  receiver  designs:   a  linear  correlation 
receiver,  the  optimum  receiver  when  the  noise  corrupting 
the  signal  is  gaussian;  a  dual  channel  design,  based 
upon  the  observed  non-gaussian  characteristics  of 
bandlirnited  low  frequency  atmospheric  noise,  which 
should  be  optimum  for  this  noise.   In  Chapter  3»  we 
presented  the  calculations  necessary  to  relate  the 
digital  simulation  to  the  analog  system  being  modelled. 
In  Chapter  4,  we  presented  our  results  for  the  two 
receivers  described  above,  and  some  bounds  upon  those 
results.   In  this  chapter,  we  will  give  our  conclusions 
based  upon  our  results,  and  we  will  also  give  some 
suggestions  for  further  research  with  the  dual  chan- 
nel receiver  design. 

5.1  Conclusions 

Our  results,  shown  in  graphical  form  in  Figures 
Jj-,1  through  '4,'4t    were  favorable.   We  found  that  the  dual 
channel  receiver  yielded  significant  performance  improve- 
ments over  the  single  channel  receiver  in  the  three  more 
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violent  weather/noise  conditions  (frontal,  tropical,  and 
quiet  night).   In  the  nearly  gaussian  quiet  noise  condition, 
the  dual  channel  receiver  performed  approximately  the 
sane  as  the  single  channel  linear  receiver*   The 
performance  improvements  in  the  more  violent  noise  cases 
were  more  than  an  order  of  magnitude  as  the  signal-to- 
noise  ratio  increased. 

The  results  were  as  we  expected.   We  note  that  the 
single  channel  linear  correlation  receiver  is  the  optimum 
gaussian  receiver,  and  it  tends  to  be  optimum  in  the 
quiet  noise  case.   However,  as  the  noise  becomes  more 
non-gaussian,  xsre  note  that  the  dual  channel  receiver 
performs  much  better.   Finally,  in  the  frontal  noise 
case,  the  most  violent  noise  case,  the  single  channel 
receiver  could  do  little  better  than  guess,  while  the 
dual  channel  receiver  made  correct  decisions  with  about 
98^  confidence  at  the  same  signal-to-noise  ratio. 

We  therefore  conclude  that  the  dual  channel  receiver, 
while  being  approximately  twice  as  complicated  as  the 
single  channel  receiver,  yielded  a  significant  improve- 
ment and  should  be  investigated  further.   It  must  be 
noted  that  we  have  simulated  only  the  simplest  case, 
using  coherent  PSK  signals,  but  we  feel  that  our  results 
indicate  that  the  dual  channel  receiver  design  is  a 
potentially  powerful  design  for  low  frequency  digital 
communication  systems. 
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5.2  Sus-p-estions  for  farther  Research 

We  have  tested  our  dual  channel  receiver  design 
against  the  optimum  receiver  for  gaussian  noise,  but 
we  have  not  tested  it  against  a  receiver  which  is  more 
similar  to  those  currently  used  to  receive  signals 
corrupted  by  non-gaussian  atmospheric  noise.   Practical 
correlation  receiver  designs  for  use  in  low  frequency 
signalling  systems  use  nonlinear  devices,  placed  before 
the  integrator,  to  limit  the  effects  of  the  non-gaussian 
high  energy  noise  bursts.   Typical  nonlinear  devices 

which  are  used,  or  suggested  for  use,  are  clippers, 

1  2 
limiters,  and  hole  punchers.    (see  Figure  5*1) 

We  suggest  that  a  possible  research  program,  based 
upon  our  results,  would  be  to  determine  the  optimum 
nonlinear  device  to  be  placed  in  the  single  channel 
receiver,  and  to  determine  the  optimum  point  at  which 
limiting,  or  attenuating.,  occurs.   Simulations  could  be 
run,  using  this  optimum  nonlinear  receiver,  and  the 
results  compared  with  the  results  we  have  obtained  with 
the  dual  channel  receiver.   Since  the  single  channel, 
nonlinear  receiver  would  be  simpler  to  implement  than 
the  dual  channel  design,  it  would  be  interesting  to 
determine  how  much  performance  improvement,  if  any ,  is 
possible  with  the  dual  channel  design  over  the  single 
channel  nonlinear  design. 

In  our  research,  we  assumed  we  had  phase  coherent 
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Figure  5»1  Nonlinear  Devices 
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receivers,  and  we  postulated  that  the  results  we  obtained 

would  give  an  indication  of  the  performance  to  be  expected 

in  the  case  of  incoherent  receivers.   Several  standard 

communication  theory  references  (e.g.  3>^)  could  be  used 

to  find  implementations  of  incoherent  receivers,  and  the 

test  we  made  could  be  redone  using  these  incoherent 

receivers. 

We  found  several  references  to  an   incoherent  receiver 

5,6 
design    which  uses  a  combination  of  energy   detection 

anc1  a  pilot  tone  to  determine  the  phase  of  the  received 

signal.   This  receiver  has  been  simulated  and  results 

have  been  presented  for  the  case  in  which  the  corrupting 

noise  is  gaussian.    It  would  be  interesting  to  implement 

this  receiver  in  a  dual  channel  design,  and  to  test  it 

in  the  presence  of  the  non-gaussian  atmospheric  noise. 

The  information  from  the  second  channel  can  also  be 

used  in  an  error  detection  scheme.   Since  the  noise  has 

a  time  varying  power  level,  the  actual  signal-to-noise 

ratio  for  each  bit  of  received  signal  may  fluctuate.   The 

second  channel  noise  envelope  measurements  can  be  used 

to  determine  the  actual  signal-to-noise  ratio  on  any  bit, 

and  when  the  actual  SNR  is  below  a  certain  threshold, 

the  receiver  can  indicate  that  the  bit  is  probably  in 

1 
error. 

Since  we  found  that  the  single  channel  receiver 

is  indeed  the  optimum  receiver  during  the  quiet  noise 
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periods,  while  the  dual  channel  receiver  performs  much 
better  in  the  more  severe  noise  periods,  it  would  be 
convenient  to  have  some  device  which  determines  when 
to  switch  between  the  single  and  dual  channel  receivers. 
Furthermore,  the  thresholds  for  the  dual  channel  design 
very  with  the  particular  weather/noise  condition  in 
existence.   Therefore,  it  would  be  beneficial  to  design 
an  algorithm  which  samples  the  noise  envelope  on  the 
second  channel  and  determines  which  receiver  and/or 
which  threshold  (v. ,  from  equation  2.16)  value  to  use. 
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